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Abstract 

The  interaction  between  electromagnetic  radiation  and  two  two-level  atoms 
with  different  resonance  frequencies  near  a  perfectly  conducting  metal  surface 
is  considered.  The  atom-atom  and  atom-surface  separations  are  assumed  to  be 
smaller  than  the  corresponding  mean  resonance  wavelength.  A  quantum-mechanical 
version  of  the  image  method  is  adopted  to  study  the  spontaneous  emission  by 
such  an  atomic  pair.  Within  the  framework  of  this  approach,  each  individual 
atom  and  its  corresponding  image  are  kinematically  correlated,  while  dynamically 
they  are  in  effect  independent.  The  total  radiation  rate  of  the  atomic  system 
is  calculated  as  a  function  of  time  for  various  values  of  the  frequency  differ¬ 
ence.  Explicit  results  are  given  for  several  different  initial  states  of  the 
atomic  system.  Some  of  them  exhibit  superradiance,  and  some  initially  act  as 
photon- trapping  states  and  eventually  are  able  to  undergo  radiative  decay. 
Oscillations  as  a  manifestation  of  beating  appear  in  the  time  evolution  of 
the  radiation  rate  in  all  cases  of  the  various  initial  states,  when  the  fre¬ 
quency  difference  becomes  larger  than  a  critical  value  given  as  twice  the 
mean  half-width  of  the  atomic  resonance  lines. 


I.  Introduction 


Over  a  decade  ago,  with  the  technical  advances  in  fatty-acid  monolayer 
1-3 

assembly,  a  series  of  measurements  was  successfully  made  on  the  fluorescent 

lifetime  of  oriented  dye  molecules  held  at  well -controlled  distances  from  a 
4-6 

metal  mirror.  It  was  found  that  for  large  distances  from  the  metal  surface 

the  fluorescent  lifetime  oscillates  as  a  function  of  distance,  while  for  small 

distances  the  lifetime  approaches  zero  monotonically. ^  A  theoretical  basis  for 

7  8 

these  experimental  observations  was  provided  by  Kuhn,  ’  who  utilized  Sommer- 

g 

feld's  classical  electromagnetic  treatment  for  radio  waves  propagating  along 
the  earth's  surface.  Within  the  framework  of  this  theory,  the  emitting  molecule 
acts  as  an  oscillating  dipole  near  a  partially  absorbing  and  partially  reflecting 
surface.  The  interference  between  the  reflected  wave  and  the  initial  wave 
gives  rise  to  the  observed  oscillations  in  the  lifetime  as  a  function  of  dis¬ 
tance.  A  number  of  researchers^”^  have  further  developed  the  calculations 
along  the  same  line  of  arguments  for  more  detailed  energy  transfers  between  the 
molecule  and  various  surfaces,  obtaining  good  agreement  with  experimental 
results.4'7-'4''5 

On  the  other  hand,  Morawitz^  adopted  a  quantum-mechanical  viewpoint  to 
investigate  the  emission  by  a  two-level  atom  at  a  distance  comparable  to  the 
radiation  wavelength  from  a  metal  surface.  Assuming  the  metal  to  be  a  perfect 
conductor,  i.e.,  with  an  infinite  conductivity,  he  replaced  the  metal  mirror 
by  an  image  behind  the  mirror,  at  a  distance  equal  to  that  between  the  source 
atom  and  the  mirror.  A  symmetric  or  antisymmetric  linear  combination  of  the 
atomic  excited  and  ground  levels  can  then  be  used  to  represent  the  initial 
state  of  the  emitting  atom  in  front  of  the  mirror,  corresponding  to,  respec¬ 
tively,  the  perpendicular  (in  phase)  or  parallel  (out  of  phase)  dipole  transi¬ 
tion  to  the  mirror  plane.  The  associated  physical  interpretation  is  simply 


that  the  emitted  photon  from  the  atom  carries  no  information  as  to  which  process 
has  occurred,  be  it  by  either  direct  emission  or  by  reflection  from  the  surface, 
so  that  it  can  be  considered  to  be  emitted  from  a  fictitious  image  atom.  The 
same  results  were  obtained  as  those  from  the  classical  approach.^® 

This  quantum-mechanical  treatment  was  later  extended  by  Milonni  and 
Knight* ^  in  the  consideration  of  spontaneous  emission  of  an  atom  between  two 

1 8 

infinite  plane  mirrors  parallel  to  each  other.  Previous  theoretical  approaches 

for  this  two-mirror  problem  usually  involve  an  expansion  of  the  electromagnetic 

field  in  appropriate  mode  functions  satisfying  the  boundary  conditions  imposed 

by  the  mirrors.  Again,  both  approaches  lead  to  the  identical  results,  lending 

the  support  to  the  quantum-mechanical  treatment  for  the  radiative  emission  of 

atoms  near  a  metal  surface.  As  a  matter  of  fact,  it  has  motivated  the  application 

20  21 

of  the  more  sophisticated  quantum-electrodynamic  theory  to  such  a  problem.  * 

In  this  paper,  we  shall  consider  the  spontaneous  emission  from  a  system  of 

two  atoms  with  different  resonance  frequencies  in  the  presence  of  a  metal  surface. 

It  is  important  to  note  that  the  spontaneous  emission  of  radiation  has  interesting 

properties  related  to  the  source  atom's  environment.  Perhaps  the  most  famous 

22 

example  of  such  an  environmental  influence  is  Dicke's  superradiance,  which  has 

stimulated  a  large  number  of  investigations,  both  theoretically  and  experiment- 

23- 29  30  31 

ally,  on  cooperative  radiation  for  atomic-gas  and  condensed-matter  * 
systems.  One  of  our  aims  here  is  to  study  a  somewhat  related  effect  of  interest, 
namely,  the  influence  of  images  on  the  radiation  of  a  system  of  atoms,  which 
may  be  in  a  stage  of  cooperative  radiation,  when  they  are  in  the  neighborhood 
of  a  metal  surface.  For  this  purpose,  the  quantum-mechanical  approach  is 
indispensable.  Another  important  feature  underlying  the  present  problem  is 
that  the  resonance  frequencies  of  the  atoms  are  different.  As  shown  in  the 
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case  of  free  space,  namely,  in  the  absence  of  the  metal,  the  emission  rate 
of  the  radiation  from  two  atoms  with  different  resonance  frequencies  has  dis¬ 
tinctively  different  patterns,  depending  on  the  ratio  of  frequency  difference 
to  the  single-particle  decay  rate.  The  situation  is  expected  to  become  more 
complicated  when  there  is  a  metal  surface  nearby.  We  shall  apply  the  image 
method^*^  to  this  problem  for  which  the  atom  and  the  corresponding  image  are 
treated  on  the  same  footing.  In  effect,  the  atoms  and  images  can  be  treated 
as  if  they  were  dynamically  independent  of  each  other,  as  far  as  the  interaction 

with  the  r  Jiation  is  concerned,  The  spontaneous  emission  rate  for  the  system 

24 

can  thus  be  calculated  by  standard  methods.  When  the  frequency  difference  is 
set  equal  to  zero,  the  case  for  two  identical  adatoms  on  a  metal  surface  ensues 
natural ly. 

This  paper  is  organized  as  follows.  In  the  next  section,  we  introduce  the 
atomic  and  radiation-field  operators.  This  is  followed  by  the  construction  of 
atom-image  correlated  states  and  the  effective  Hamiltonian  for  the  system,  where 
physical  arguments  are  provided.  The  equations  of  motion  are  given  in  Sec.  III. 
Here  a  conservation  law  is  established  which  is  useful  for  expressing  the  total 
emission  rate  solely  in  terms  of  the  atomic  operators,  and  subsequently  a 
hierarchy  of  coupled  equations  are  solved  with  help  of  appropriate  approximations. 
Finally,  we  are  able  to  relate  the  emission  rate  to  all  expectation  values  of 
initial  states.  For  Sec.  IV,  a  variety  of  possible  initial  states  which  can  be 
prepared  experimentally  are  considered,  and  the  corresponding  emission  rates 
are  computed  with  respect  to  time  and  frequency  differences,  accompanied  by  a 
discussion  of  the  results. 

II.  Theory 

A.  Two-Level  Atoms  and  the  Quantized  Electromagnetic  Field 

We  shall  restrict  ourselves  to  the  case  where  each  atom  near  a  metal  surface 
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has  just  two  energy  levels.  We  let  1-  >  denote  the  ground  level  for  the  atom 
labeled  by  the  index  p  and  l+y>  the  excited  level,  where  normalization  and 
othogonality  require  that<+y|+y>  =  <-y|-  >  =  1  and<+y|-y>  =  ^e  designate 
the  frequency  corresponding  to  the  energy  difference  between  the  two  levels 
is  u  ,  where  Although  we  shall  concentrate  on  the  situation  in  which 

the  distance  between  the  atoms  is  much  smaller  than  their  radiation  wavelengths, 
we  shall  assume  that  the  wave  functions  of  the  two  atoms  do  not  overlap,  i.e., 
<+y|±v>  =  0  for  pf  v,  so  that  the  atoms  are  distinguishable.  The  transitions 
between  the  ground  and  excited  levels  can  be  facilitated  by  introducing  excita- 

4. 

tion  (raising)  and  de-excitation  (lowering)  operators  cy  and  cy,  respectively  , 
with  the  properties33  that  cjj-u>  *  |+y>,  cy|+y>  =  0,  cy|-y>=0  and  cy|+y>  =  |-y>. 
As  a  result,  it  can  be  readily  recognized  that  cy  and  cy  obey  the  anti  commutation 
relation,  i.e.. 


{  VcJ)a  1-  (2.1 ) 

It  is  also  required  by  the  nonoverlapping  property  of  the  two  atoms  that  operators 
representing  different  atoms  commute,  i.e.. 


’Cv]  =  CVCv]  =  °’ 


(2.2) 


[Cy’Cy]"  O’  f0r  v  *  V' 

Thus  the  atoms,  according  to  the  above  descriptions,  can  be  modelled  as  dipoles. 

Since  our  main  concern  here  is  the  total  radiation  rate  of  the  system,  we  shall 

neglect  the  dipole-dipole  interaction  between  the  atoms,  which  is  expected  to 

34 

only  shift  the  emitted  radiation  frequencies.  Therefore,  the  atoms  are  inter¬ 
acting  with  each  other  via  the  common  radiation  fields. 

Each  mode  of  the  radiation  field  can  be  treated  as  a  quantum  oscillator  with 
frequency  fl  .  We  further  distinguish  the  n^  different  quantum  oscillators  with 
the  same  frequency  by  a  second  index  8.  For  example,  in  the  case  of  a  plane 
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wave,  a  may  refer  to  the  frequency  (or  |ic|,  where  £  is  the  wave  vector),  while 
6  refers  to  direction  of  £  and  the  polarization.  Transitions  between  occupation 
number  states  |nag>  are  described  by  annihilation  and  creation  operators  a^  and 

aT„  with  the  well-known  properties,  a  Jn  _>  =  /n  Jn  c-l>,  and 

a6  aB1  aB  aB1  aB 

aotB I na.6>  =  !nr,ft+l>-  These  operators  obey  the  commutation  relations 


aB  '"aB 


^aaB’aa'B'^  =  6aa'6BB1’ 


(2.3) 


and  the  occupation  number  states  |nag>  form  the  orthonormal  set,  namely, 

<n  Jn  ,„i>  =  6  All  atomic  operators  c  and  c+  commute  with  all  operators 

aB1  a  B  aa  B6  u  u 

(a^g  and  a^g)  for  the  quantized  electromagnetic  oscillator. 

B.  Atom-Imaqe  Correlation  in  the  Presence  of  a  Metal  Surface 


In  the  present  work,  we  shall  assume  the  metal  to  have  a  perfect  conductivity. 
Therefore,  the  role  played  by  the  metal  surface  with  respect  to  radiation  is  nothing 
more  than  a  reflecting  mirror.  The  radiation  emitted  from  an  atom  close  to  the 
surface  can  reach  another  atom  by  either  direct  transmission  or  through  reflec¬ 
tion  from  the  surface,  as  though  it  were  emitted  by  the  image  of  the  first  atom 

(see  Fig.  1).  There  is  therefore  a  possibility  of  interference  between  the 
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radiation  fields  from  two  processes,  as  in  the  case  of  a  Lloyd's  mirror.  Since 
these  fields  carry  no  information  about  which  process  has  occurred,  it  leads 
one  to  consider  the  atom  and  the  corresponding  image  on  the  same  footing.  7 
Namely,  we  may  describe  the  emission  by  a  single  atom  in  the  vicinity  of  a  metal 
surface  as  a  two-atom  problem  with  complete  uncertainty  as  to  which  atom  is 
excited.  We  can  thus  write  the  initial  (excited)  state  incorporating  such 
ambiguity  for  atom  p  as 


ivt  ■  -y  !■ vhw*1' v'V' 


(2.4) 


where  the  indices  A  and  I  refer  to  the  atom  and  image,  respectively,  while  the 


meanings  of  l+>  l->  are  the  same  as  those  given  in  Sec.  II. A.  State  |e  >+  is 
used  if  the  atomic  transition  dipole  is  perpendicular  to  the  plane  of  the  metal 
surface,  because  the  dipoles  of  the  atom  and  its  image  have  the  same  phase.  On 
the  other  hand,  state  Je^>_  is  used  when  the  transition  dipole  is  parallel  to 
the  plane  of  the  metal  surface  where  a  dipole  has  a  phase  opposite  to  that  of 
its  image.  The  physical  picture  embedded  in  Eq. (2 .4)  represents  a  quantum  mechan 
ical  version  of  the  classical  description  of  an  oscillating  dipole  and  its  image. 
The  ground  state  is  given  by 


IV  =  IV'V- 


(2.5) 


C.  Hamiltonian  and  Emission  Rate 


From  the  preceding  discussion,  we  see  that  the  effect  of  the  metal  surface 
can  be  cast  in  terms  of  the  correlation  between  the  atoms  and  their  corresponding 
images.  The  latter  can  be  treated  effectively  as  independent  but  identical  to 
their  respective  source  atoms,  in  the  sense  that  the  source  atom  and  its  image 
are  kinematically  correlated  but  in  effect  dynamically  independent.  Thus,  one 
can  immediately  write  down  an  effective  Hamiltonian  for  the  system  of  two  atoms 
with  different  frequencies  and  near  a  metal  surface  as 


2  t  t 

H  r-  =  11  E  w  ( c,  c.  +CT  c,  ) 
eff  =|  p'  Au  Au  Iu  Iir 


n 


-  HgE  +  (cA+u+ci>0(J}, 

M  CL  p 


(2.6) 


where  g  is  the  radiation-atom  coupling  constant  and  c.  (c,  )  is  the  creation 


(annihilation)  operator  for  p-th  atom  and  cju(cjy)  for  the  corresponding  image, 
which  obey  the  operator  algebras  given  in  Sec.  II. A.  In  writing  Eq.(2.6)  we 
have  assumed  that  the  atom-atom  and  atom-surface  separations  are  both  smaller 
than  any  of  the  radiation  wavelengths,  so  that  there  is  no  phase  difference  at 
different  positions  of  atoms  and  images.  Eq . (2.6)  in  fact  describes  effectively 
the  same  physical  properties  as  those  of  four  atoms,  two  of  them  with  resonance 
frequency  and  the  other  two  with  a  different  resonance  frequency  u^,  except 
that  two  with  the  same  frequency  are  correlated  to  each  other  rather  than  com¬ 
pletely  free  atoms. 

The  quantity  of  central  interest  is  the  total  emission  rate  of  the  system, 
r(t),  which  is  given  by 

r(t>  -It“<al8(t)ac,8!t»> 

a  p 

where  a1”  and  a  ,  are  expressed  in  the  Heisenberg  picture.  [From  hereon  for 

CXp  up 

simplicity  of  presentation,  the  time  dependence  of  the  operators  will  not  be 
displayed  unless  necessary.]  The  advantage  of  working  in  the  Heisenberg  picture 
is  that  the  expectation  value  in  Eq.(2.7)  is  taken  on  the  initial  state  of  the 
whole  system,  |c(0)>, which  consists  of  the  atomic  and  radiation-field  states. 
Since  we  are  only  interested  in  the  spontaneous  emission,  it  is  understood  that 
the  radiation-field  part  of  the  state  involved  here  is  the  vacuum  state.  Thus, 
we  only  have  to  pay  attention  to  the  various  atomi c  initial  states  of  the  system. 

III.  Calculations 
A.  Equations  of  Motion 

Using  the  effective  Hamiltonian  of  Eq . (2.6)  and  Heisenberg's  equation  of 


motion,  one  can  obtain  the  following  equations  of  motion: 
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”Au 


iu>  c. 

U  Ay 


=  i g  [c.  ,ct  ]E  Z  a  0 , 

L  Au  Ay J  Q  aB 
a  B 


(3.2) 


'Iy 


lai  CT 

y  Iy 


1  9  £cIy,cIy 


Ja  Z  3  , 

a  B  aB 


(3.3) 


where  the  dot  over  an  operator  signifies  the  time  derivative.  In  order  to  simplify 
the  calculation,  we  begin  by  eliminating  the  rapidly  oscillating  part  in  the  opera¬ 
tors  by  choosing  a  frequency  o)q  between  u>-j  and  u)^  and  defining  the  operators 
A  and  C  through  the  relations 


;a  fl(t)  =  A  (t)e  0 
aB  a 

(3.4) 

(3.5) 

-i«*>nt 

yt)  -  C,u(t)e  V 

(3.6) 

These  new  operators  are  clearly  less  oscillatory  in  time.  To  write  Eqs.(3.1)  - 
(3.3)  in  terms  of  these  new  operators,  we  sum  (3.1)  over  B, 

Za+iflia=ign  l(c.  +cT  ) ,  (3.7) 

0  aB  a  aB  a  a  '  Ay  Iy' 

B  a  y 

and  then  transform  it  by  using  Eqs . ( 3 . 4 ) - ( 3.6)  into 

\  +  1(V"0)Ac,  =  lgnc.J  'W1  (3'8) 

where  we  have  used  the  fact  that  there  are  n  modes  of  radiation  oscillators 

a 

having  the  same  frequency  Similarly,  we  have  from  Eq . (3.2) 

%  *  Aa  <3'9> 

a 


and  from  Eq . (3. 3) 


Cly  +  *'Wy_<tyCIy  =  igLCIy’CIy^  V 


(3.10) 


Eq . ( 3. 8)  can  be  integrated  to  give 


-i(fl -o)n) 

A  =  e  a  0  A  (0) 
a  a'  ' 

t  1(w  -ui0)  ( t*-t) 

+  i9n“  ^dt'e  CVCIu3 


(3.11) 


Summing  over  all  a,  we  obtain 


where 


-i(fl  -wn)t 

H  =  I  e  01  0  A  (0)  +  igF(t), 
a  a 


t  i(fl  -wn)(t'-t) 

F(t)  =  p/0dt,CVCIu]^ap(Qa)e 


(3.12) 


(3.13) 


and  p(wa)  is  the  density  of  radiation  oscillators  per  unit  frequency.  Following 

24 

the  arguments  of  Dillard  and  Robl  that  the  interaction  is  confined  to  some 
frequency  range  wq  -  <_  Sl<_  uig  +6w  with  6u>  «  Wg,  the  double  integrals  in 

Eq .(3.13)  can  be  performed  to  yield 


F(t)  =  TTP(a,0)E  (CAy+CIy), 


(3.14) 


Consequently,  Eq . (3.12)  becomes 


I  A  =  z  e 

a 

a  a 


Aa(0)  +  iTrgp(u)0)I  (CA  +Cj  ) 
U 


(3.15) 


or,  with  Eqs .(3.4) -  (3.6), 


“  E/aS(0)  + 

a  6  a  6  y 


(3.16) 


At  this  stage  we  are  ready  to  obtain  the  total  emission  rate  by  determining 
a  means  of  finding  the  expectation  value  of  a^a^  [see  the  right-hand  side  of 
Eq . (2.7)3.  We  multiply  Eq .(3.1)  on  the  left  by  atQ  and  the  Hermitian  adjoint  of 


Eq. (3.1)  on  the  right  by  a^,  and  the  two  results  are  summed  to  give 

it  “  aIsaae  =  19a  eJ{aQ‘6(cAu+CIli)  '  (cAu+CIu)aaB} *  (3>17) 

Substituting  Eq.  (3.16)  and  its  Hermitian  adjoint  into  Eq .(3.17)  and  taking  the 
expectation  value  on  the  initial  state,  we  obtain 

■k  “IbW  =  2,192  pt“0  )lZ«clv  +cJ’„)(cAv+cIv>  ’ 
a  8  jj  v 

'  V^<(cAm+CIu)(cAv+cIv)>-  (3.18) 

p 

where  y  =  2 ng  p(wg)  is  the  transition  rate  for  spontaneous  emission  by  an  isolated 
excited  atom.24’37  Eq .(3.18)  can  be  further  wri tten  as 

dt  Jg<aIeaaB>  =  y^<cAu%  +  CIy  CIy  > 

+  YJ^<{i+ClJ{cAv+CIv)> 

+  Y^CA1CI1+CA2CI2+CI1CA1+CI2CA2^ ’  (3,19) 

whose  solution  then  yields  the  total  emission  rate. 


B.  Conservation  Law 


We  digress  for  a  moment  here  to  derive  a  conservation  law,  which  will  prove 
to  be  helpful  solving  the  overall  problem.  From  Eq.(3.2)  and  its  adjoint,  one 
can  write 


dt  I  CAyCAy  =  192  CAytcAy’cAy^  aa6  "  1  9  Z  l  Sa6  ^cAy’cAy^CAy  ’  {3'20> 

U  y  aB  a  8  y 


From  Sec.  II. A,  we  know  that  cJy[cAvj,cJu]  =  cAu  and  c^u.  As  a 


consequence,  Eq.(3.20)  becomes 


(3.21] 


(3.22] 


-jT-ici  c»  =  igZc:zzao-igZZa0Zc.. 
dty  Au%  %A,a6a 8  aa6  «6y  Au 

Similarly,  from  Eq.(3.3)  and  its  adjoint,  we  obtain 


J  j  x  J. 

57  Z  cT  cT  =  1  g  Z  cT  Ila  -  igZZa  fiIc.  . 

dt  y  IU  IU  P  lMa6  aB  <xB  a6U  IU 


The  addition  of  Eqs.(3.21)  and  (3.22)  gives 


?t  l(cAyCAy+CIyCIu)  =  1  9  11  l(cAy+CIy)aaS  *  1  9  H  [  *l&{cbv+Cl\?  *  (3‘23 

p  Ot  p  U  u  p  p 

In  comparison  with  Eq . ( 3- 1 7) ,  we  find  the  following  conservation  law: 


Jt<VVNI>  •  °- 

Nph  rePreSents  the  total  number  of  photons,  defined  as 


N.  =  I  1  d  Q  d  q  $ 
ph  a&  uB  aS 

represents  the  total  number  of  excited  atoms,  defined  as 


NA  '  ICAyCAy’ 

u 


and  represents  the  total  number  of  excited  images,  defined  as 


N,  -  Zc,  c.  . 

I  u  Iu  In 


(3.24 


(3.25 


(3.26 


(3.27 


Eq.(3.24)  states  that  the  total  number  of  the  photons  and  the  excited  atoms  and 
images  is  constant.  Note  that  it  is  an  operator  equation  and  therefore  holds 
for  expectation  values  on  any  states.  In  other  words,  Nph  +  +  Nj  is  a 

constant  of  motion. 


r . 


C.  Solutions 

Having  set  up  the  conservation  law,  we  immediately  recognize  that,  from 
Eqs . (3.24)  -  (3.27), 


■jr  E  z  a+Qa  Q 
dt  a8  aB  aB 


dt  ^cAy%+cIucIy^ 


(3.28) 


Thus,  with  Eq.(2.7)»  we  have  the  total  radiation  rate  given  by 


r(t)  *=  -  Z  <camcAm+cImCIm>- 


(3.29) 


Furthermore,  by  combining  Eqs. (3. 19)  and  (3.28)  we  have  an  equation  involving 
only  the  atomic  operators. 


-rr  Z  <  c|  c.  +ct  CT  >  +  yZ  <ct  C.  +ct  CT  > 
dt  Am  Am  Im  Im  t  Am  Am  Im  1m 


=  -^<(CWm)(CAv+CiJ>  "  Y<CA1CI1+CA2CI2+CI1CA1+CI2CA2>  (3,30} 

At  this  point,  it  should  be  emphasized  that  the  expectation  value  of  off- 
diagonal  terms  on  the  right-hand  side  of  the  above  equation  does  not  necessarily 
vanish.  This  is  the  place  where  the  coherence  comes  into  play.  It  is  in  the 
form  of  a  dipole-dipole  interaction  between  atoms  and  images,  as  induced  by  the 
common  radiation  with  which  the  atoms  interact  in  phase,  resulting  from  the  fact 
that  the  atom-atom  and  atom-surface  separations  are  smaller  than  the  mean  resonance 
wavelength.  Were  we  to  include  the  Jirect  dipole-dipole  interaction,  we  would 
have  an  additional  term  of  the  same  form  but  with  a  direct  dipole-dipole  coupling 
coefficient  other  than  y  on  the  right-hand  side  of  Eq.(3.30).  However,  one  should 
note  that  the  expectation  value  of  the  direct  dipole-dipole  interaction  vanishes, 
since  it  only  has  nonvanishing  off-diagonal  matrix  elements  due  to  the  selection 
rule  of  electric  dipole  transitions.  The  direct  dipole-dipole  interaction  is 
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therefore  expected  to  affect  the  frequency  shift  only. 

To  solve  Eq.(3.30),  we  define  M.  P  and  Q: 

M  =  jj(cAucAy+cIucIyJ  (3,31) 

P  =  CA1 CA2+CA1 CI2+CI 1 CA2+CI 1 CI 2  *3‘32) 

and 

Q  =  CA1CI1+CA2CI2+CI1CA1+CI2CA2 ‘  (3-33) 

Thus,  Eq .(3.30)  becomes 

m  +  ym  =  -  y(P+P  +q)  (3.34) 

where  the  letters  in  lower  case  represent  the  expectation  values  of  the  corres¬ 
ponding  operators  of  capital  letters  defined  in  Eqs .  (3.31 )  -  (3.33),  i.e., 
p  =  <  P>,  etc. 

In  order  to  obtain  p  and  q,  we  consider  the  equations  of  motion  satisfied 
by  P  and  Q, 

P  +  (2y+iA)P  =  -y(M  +  Q  -  cJ1cA1(cJ2+c}1+c|2)(cA2+cI2) 

"  CI1CI1(CA1+CA2+CI2)(CA2+CI2) 

^CA1+CI1  ^CA1+CA2+CI1  ^CI2CI2 

^CA1+CI1^CA1+CI1+CI2^CA2CA2}  ^3'35) 

and 

q  +  q  =  -Y{M+P+P+-cJ1cA1(c^2+c|1+c}2)cI1 

■  CA2CA2^CA1+CI2+CI1 ^CA2  "  CI1CI1 (cAl+cA2+cI2 )cAl 
-  c{2cI2(cA1+cA2+c{1)cI2  -  cAi(cAl+cA2+cI2^cIlCn 

"  CA2^CA1+CA2+CI1)CI2CI2  ~  CI1 ^CA2+CI1+CI2)CA1CA1 
’  CI2^CA1+CI1+CI2)CA2CA2}* 


(3.36) 


where  A  *  ^  •  uj . 

The  off-diagonal  terms  of  four  multiplicative  atomic  operators  in  the  last 

•  +  + 

two  equations,  such  as  CA2CI2  etc-»  are  non-vanishing  provided  that  the 
two  atoms  exchange  a  photon  with  a  third  one  in  its  excited  level.  If  we  neglect 
such  terms,  it  is  equivalent  to  saying  that  we  are  neglecting  small  terms  on  the 
order  of  -*-<  1.  Under  this  approximation  we  have 


and 


P  +  (2y+iA)P  =  ->{M+Q-2R1} 


Q  +  YQ  =  -Y{M+P+P-2R2), 


(3.37) 

(3.38) 


where 


R1  =  (cAlCAl+cIl  CI1^CA2CA2+CI2CI2)  (3l39j 

and 

R2  =  2(cAlcAlCnCn+CA2CA2cI2CI2)‘  (3'40> 

Here  we  have  made  use  of  the  fact  that  CcAacAa,cIcxCIc^  =  0  etc’  t0  write  Ri  and 
R2 ,  respectively,  in  the  compact  form  in  Eqs.(3.39)  and  (3.40).  The  corresponding 
equations  for  the  expectation  values  of  Eqs .(3.37)  and  (3.38)  are 


p*  +  (2y+i'A)p  =  -Y(m+q-2r1)  (3.41) 

and 

q  +  Yq  =  -Y(m+p+p*-2r2).  (3.42) 


From  Eq.(3.41),  we  further  obtain 

s  +  2ys  +  iAd  =  -Y(2m+q+q*-2r1-2r^*) 


(3.43) 


and 


d  +  2yd  +  iAs  =  -2y[ri  -  r,) 


(3.44) 


where  s  =  p  +  p  and  d  =  p  -  p*.  Taking  the  time  derivative  of  Eq . (3.43) ,  we  get 


s  +  2ys  +  iA3  =  -Y(2m+q+q*-2r-|-2r-j*) .  (3.45) 

We  note  here  that  q=q*,  since  q  and  q*  satisfy  the  same  equation  (3.42),  due  to 
the  fact  that  Rg,  as  defined  in  Eq. (3.40) ,  is  a  Hermitian  operator  and  therefore 
r2=  r^*.  Thus  we  can  rewrite  Eq.(3.45)  as 

s  +  2ys  +  i  =  -2Y(m+q-r-]-r^*) .  (3.46) 

Substituting  Eqs.(3.34),  (3.42)  and  (3.44)  into  Eq.(3.46),  we  have 

s  +  2ys  +  A2s  =  2y{iAd  +  iA(r^*-r-j)  +  2y(m+s+q)  -  2Yr2+^i+^|*  ^  •  (3.4 

Furthermore,  using  Eqs.(3.43)  and  (3.47)  to  eliminate  iAd,  we  obtain 

s  +  4ys  +  A2s  =  2y{  2y(r1+r1*)+r-|+r1*+iA(r1*-r1  )-2yr2  }.  (3.48) 

As  shown  in  the  Appendix,  r^  =  r-|(0)e"2Yt  and  r2  =  r2(0)e~2Yt.  Equation 
(3.48)  then  becomes 

s  +  4ys  +  A2s  =  -4y2r2(0)e  2y*.  (3.49) 

To  simplify  the  above  expression,  we  further  introduce  the  dimensionless  variables 

p  2  2 

<t>=yt  and  A  =  A  /y  .  Eq . (3.49)  then  takes  the  new  form  as 

s"U)  +  4s'  (<{>)  +  A2s(4>)  =  -4r2(0)e’2<t>,  (3.50) 

where  the  primes  refer  to  the  derivatives  with  respect  to  <j>.  The  general  solution 
of  this  equation  is  immediately  given  by 

-£+4>  -C.41  4  _oA 

s($)  =  K,e  +  K.e  + - T  r„(0)e  (3.51) 

1  4 -A*  L 

where  i  -  2  *  except  for  the  special  case  A  =  2  which  gives  rise  to  a 
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divergent  solution.  and  are  determined  by  the  initial  conditions  s(0) 
and  s'(0)  and  will  be  given  later  on  for  the  different  cases  of  X. 

Having  obtained  the  solution  for  s(<ji),  we  now  proceed  to  solve  the  other 
equations.  In  terms  of  the  new  variable  4>»  Eqs . (3.29) ,  (3.34)  and  (3.38)  become, 
respecti vely, 

r(t)  =  -  ym'U),  (3.52) 

n'U)  +  m(4>)  =  -  (sU)  +  q(4>)}  (3.53) 

and 

q'(<j>)  +  q(<|>)  =  ~{m(4>)  +  s(<j>)  ^r^Oe-2*} .  (3.54) 

Our  goal  is  to  find  out  m'(<f>).  For  this  purpose,  we  subtract  Eq.(3.54)  from 
Eq.(3.53)  and  obtain 

^On(<j>)  -  q(0}}  =  -2r2(0)e"24>,  (3.55) 

which  gives  the  result 

q(<t>)  =  m(4.)-m(0)  +  q(0)  +  (l-e_2<V2(0)  •  (3.56) 

Substituting  Eq.(3.56)  into  Eq.(3.53),  we  have 

m'U)  +  2m(4>)  =  -{s(<^)-m(0)+q(0)+(l-e'2,1,)r2(0)  } .  (3.57) 

Thus,  the  solution  of  Eq.(3.57)  can  be  obtained  straightforwardly. 

The  properties  of  s(<}>)  as  expressed  in  Eq .(3.51)  depend  on  the  values  of 
the  parameter  X.  Let  us  consider  the  following  cases  (i),  (ii)  and  (iii): 


(i)  0  £  X  <  2 
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2  1/2 

Writing  4  =  (4-A  )  ,  we  first  determine  and  from  the  initial  con¬ 

ditions  s(0)  and  s ' (0) : 

Ki = "  h{  ?  r2(0)  +  (2-^s(°)+  s'(°)> 

and 

K2  ^  {\  r2(0)  "  (2+^s(°)  -  s '  (0)} » 
where  s'(0)  can  be  obtained  directly  from  Eq . ( 3.43) 

s '  (0 )  =  -2m(0)  -  2q(0)  +  4^(0)  -  2s (0)  -  i  (4-42)1/2d(0) . 

Thus,  Eq . (3.57)  can  be  integrated  with  the  expression  (3.51)  given  for  A(<j>) 
and  Eqs.(3.58)-  (3.60).  The  desired  quantity  can  then  be  obtained  as 

-m'(d»)  =  e‘%^)m(0)  +  (^-)q(O)  +  \  r,(0)  -  d(0)] 

4  4  4  V 

2 

+  [4jr  ~  2 ( ~ ) <t>3 rp ( 0 )  +  [-^m(0)  +  4><l(0)  -  r^(0) 

4  4  4  4  4 

2  1/2 

-  4-  r?(0)  +  s(0)  +  i-^-i- —  d(0) ]cosh  4$  - 

4  4 

[|  m(0)  +  |  q(0)  -  |  r-j (0 )  -  ^3  r2(0)  +  2  s(0) 

2  1/2 

+  j  J - d(0) ]sinh4<j>) .  (3.61) 

(ii)  A  >  2 

Writing  n  =  (A2-4)1//2,  we  can  use  the  above  results  directly  with  4  replaced 


(3.58) 

(3.59) 

(3.60) 


by  in: 


-m' 


U)  =  e‘2*{3-^  m(0)  t  (3-ii)  q(0)  -  ^  r,(0) 


-  i  0(0)]  ♦  [-  4  -  2(^)*]r, 


(0) 


+  t-  4  m(0)  -  \  q(0)  +  4  r,(0)  +  4  r  (0)  ♦  s(0) 

n  n  n  n 


2  1/2 

1  ^4 - d(0)]cosn*-  Cfjm(O)  +  \  q(0)  -  4  r,(0) 


2  1/2 

*  4  r2(0)  +  f  s(0)  •  <  n1  d(0)3s1nn<t.)  . 


(3.62) 


(iii )  X=  2 

For  this  special  case,  we  have  to  reconsider  the  solution  s(4>)  for  Eq.(3.50), 
where  the  characteristic  polynomial  has  a  root  of  double  multiplicity  and  the 
inhomogeneous  term  has  the  same  exponential  of  the  homogeneous  solution  as  well. 


Accordingly,  the  general  solution  is  sought  to  be 
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s(<t>)  =  (K3+K44>)e'2*-  2r2(0)4>2e'2<t>, 


(3.63) 


where  and  are  again  to  be  determined  by  the  -initial  conditions  s(0)  and 
s'(0).  Following  the  same  procedures  as  before,  we  have 


-m'U)  =  e'2<^{m(0)  +  q(0)  +  s(0)  -  2[m(0)  +  q(0)  +  s ( 0 )  -2r-j(0) 


-  r2(0)  +  id(O)]0  +  2[m(0)  +  q(0)  -2^(0)  -  r2(0) 

+  id(0)]4>2  +  f  r2(0)4>3}.  (3.64) 

The  total  radiation  rate  r(t)  for  the  system  is,  as  given  by  Eq.(3.52), 
proportional  to  the  quantity  -m'(4>)»  which  is  seen  for  every  case  to  be  deter¬ 
mined  by  the  initial  values  m(0),  q(0),  r^(0),  r2(0),  s (0)  and  d(0).  These 
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quantities,  being  the  expectation  values  of  atomic  operators  on  the  initial  state 
of  the  whole  system,  depend  on  the  initial  states  of  the  system  and  are  considered 
for  a  variety  of  different  initial  states  in  the  next  section. 

IV.  Numerical  Results  and  Discussion 

The  following  initial  states  for  the  system  in  which  both  atoms  have  transition 
dipoles  oriented  perpendicular  to  the  metal  surface  are  chosen  for  specific  consideration 

(a)  |H'(0)>=  |e1>|e2>,  (4.1) 

where  both  atoms  are  initially  in  their  excited  states,  and 

(b)  |4'(o)>  =  I e1  > | g2> *  (4.2) 

where  one  atom  is  in  its  excited  state  while  the  other  is  in  its  ground  state. 

Without  loss  of  generality,  we  focus  on  the  case  in  which  the  atom  with  reso¬ 
nance  frequency  is  in  its  excited  state  and  the  one  with  is  in  its  ground 

state.  The  solution  will  be  identical  to  |4'(0)>  =  |g-|>|e-|>  since,  as  can  be  seen 

o 

from  Eqs.(3.6l),  (3.62)  and  (3.64),  it  depends  only  on  A  . 

(c)  | VF( 0 ) >  =  —{  | e1  > | g^>  +  | e2> 1 9-j >> »  (4-3) 

where  only  one  atom  is  in  its  excited  state  while  the  whole  system  is  in  a  sym¬ 
metric  combination. 

(d)  | V(0)>  =  —{  |e]>|g2>  -  |e2>|g]>}  (4.4) 

where  only  one  atom  is  in  its  excited  state  while  the  whole  system  is  in  an 
antisymmetric  combination. 

(e)  |4»(0)>  =  [4:  (|ei>  +  |gi>)]*  [—  (|e->  +  |g?>)3,  (4.5) 

/2  &  c 

where  each  atom  of  the  system  is  in  an  in-phase  superposition  of  states. 
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i 

E 


(f)  |4'(0)>  =  [—  ( I  e-.  >+ 1  g]  >)  ]  x  [^-  (|e2>-|g2>)] 
ft  ft  c  c 


(4.6) 


where  the  atom  with  frequency  is  in  an  in-phase  superposition  of  the  excited 
and  ground  levels  while  the  one  with  frequency  is  in  an  out-of-phase  super¬ 
position  of  levels.  For  the  same  reason  as  in  case  (b),  the  solution  is  identical 

to  the  initial  state,  |4'(0)>  =  [—  ( | e-. >- 1 g-. >) ]  *[—  ( |e9>+|g7>)]. 

ft  '  '  ft  c  c 

i6 


(g)  1^(0)  >  =  ^  (|g>+e  y|eu>). 


v  =  1,2 


(4.7) 


where  each  atom  is  in  an  equal  and  random-phase  superposition  of  levels.  The 
density  operator  representing  this  situation  is  thus  given  by 

P(0)  =  — -p  /d61/2d6J|’MQ)x4'1(0)|]  x  [ \T  (0)><  S'  (0)  |  ] 

00^  *  c 

=  Ljlg1><g1  l+  ^|e1><e1  |]x  [l|g2><g2|  +  i|e2><  e2|]  (4.8) 


In  this  case,  the  expectation  value  of  an  operator  o  is  interpreted  as 


<  0>  =  Tr(p(0)o).  (4.9) 

We  list  in  Table  I  all  the  expectation  values  necessary  for  evaluating  the 

emission  rate,  r(t).  The  normalized  emission  rate,  r(t)/y,  for  the  different 

initial  states  of  Eqs .(4.1 )  -  (4.7),  is  plotted  in  Figs.  2-8  with  respect  to  the 

normalized  time,  <p  =  yt,  and  the  parameter  characterizing  the  frequency  difference, 

X  =  A/y.  The  special  case  of  two  identical  atoms  is  represented  by  the  curves 

with  X=0.  These  results  can  be  roughly  divided  into  two  regions  of  different 

behavior  for  the  time  evolution  of  the  emission  rate.  For  A<2y  the  emission 

rate  exponentially  decays  towards  zero,  while  for  A  >  2y  the  rate  shows 

oscillations  before  decaying  to  zero.  We  see  that  the  dividing  frequency  occurs  at 

A  =  2y.  The  physical  explanation  can  be  given  as  follows.  For  an  excited  atom 

with  its  transition  dipole  perpendicular  to  the  metal  surface,  the  half-width  is 
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doubly  broadened  as  in  the  case  of  Dicke's  superradiance,  since  it  is  described 


by  the  state  |e>+  of  Eq.  (2.4)  as  a  result  of  phase  coherence  between  the  atom  and 
its  image.  Therefore,  if  A  <  2y,  the  two  atoms  are  hardly  distinguishable  from 
each  other,  and  they  undergo  radiative  emission  as  if  there  were  two  identical 
atoms.  Only  when  A  >  2y  does  the  difference  of  the  two  resonance  frequencies 
become  pronounced, and  beating  takes  place  as  manifested  in  the  oscillation  of  r(t). 

Let  us  now  examine  closely  the  initial  values  of  the  emission  rate, 
r(0),  which  are  used  as  a  measure  of  the  radiative  intensity  in  the  usual  sense, 
in  order  to  gain  a  better  understanding  of  the  problem.  At  t=0,  F(0)  can  be 
immediately  written  from  Eqs.(3.29)-  (3.34)  in  the  form 

r(0)  =  yi  «Jm(0)cAii(0)^m(0)cIu(0)>  ♦  yj<c;u(D)+cIu(0)][cAv(0)+cIv(0)]> 

+  Y<cJ1(0)cI1(0)+c^2(0)cI2(0)+cI1(0)cJ1(0)+c|2(0)cA2(0)> 

=  y[m(0)+q(0)+s(0)].  (4.10) 

Clearly  this  equation  contains  only  the  initial  excitation  stored  in  the  atomic 
system  described  by  m(0),  and  the  induced  dipole-dipole  interaction  among  the  atoms 
and  images  due  to  the  energy  fluctuation  of  the  vacuum  as  specified  by  q ( 0 )  and 
s(0).  Therefore  q(0)  and  s(0)  together  signify  a  certain  degree  of  coherence  in 
cooperative  emission.  It  is,  however,  important  to  distinguish  the  difference  of 
physical  natures  embedded  in  the  correlations  expressed  by  q(0)  and  s(0).  q(0), 

as  seen  through  the  definition  of  its  corresponding  operator  Q  given  in  Eq . ( 3 . 33) , 
represents  the  correlation  between  an  individual  atom  and  its  own  image.  There¬ 
fore,  q(0)  is  solely  due  to  the  existence  of  the  metal  surface.  In  contrast,  s ( 0 ) 
is  attributed  completely  to  the  correlation  brought  about  by  the  radiation  fields. 

As  time  evolves,  correlations  of  higher  orders  are  literally  needed  to  determine 
the  emission  rate.  As  we  have  seen  in  Sec.  III.C,  the  higher  correlations  are 
manifested  in  the  appearance  of  operators  such  as  R-j  and  R  ,  in  which  more  atomic 


operators  are  clustered  together. 

From  Table  I,  we  see  that  only  the  initial  state  (a)  has  2  units  of  excitation 
stored  in  the  atomic  system,  i.e.,  m{0)=2,  since  this  state  has  both  atoms  in 
their  excited  states,  while  in  each  of  the  other  cases,  (b)  -  (g),  only  one  atom  is 
in  effect  in  an  excited  state.  Similarly,  since  the  transition  dipole  can  induce 
an  image  dipole,  only  the  initial  state  (a)  which  contains  two  excited  atoms  has 
q(0)=2,  but  in  cases  (b)-  (g)  each  has  the  value  q(0)=l.  We  see  that  the 
enhancement  in  state  (a)  is  solely  due  to  the  existence  of  the  metal  surface, 
because  s(0)  =  0  as  shown  in  Table  I.  For  the  initial  state  (b),  in  which  one 
atom  is  excited  while  the  other  one  in  its  ground  state,  there  is  an  enhancement  in 
the  emission  rate  due  to  the  existence  of  the  metal  surface  as  well,  but  with  a 
reduced  factor  due  to  the  fact  that  only  the  excited  atom  can  induce  an  image 
to  participate  in  the  radiative  process. 

For  the  initial  state  (c),  in  which  the  two  atoms  are  in  an  in-phase  combina- 
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tion  resembling  the  Dicke's  superradiati ve  state,  we  have  the  highest  enhancement 

factor  in  s (0 ) ,  and  combined  effect  of  q(0)  and  s ( 0 )  results  in  an  initial  emission 

rate  r( 0)  as  high  as  in  the  case  of  the  initial  state  (a).  In  contrast,  the  initial 

state  (d),  in  which  the  two  atoms  are  in  an  out-of-phase  combination,  gives  r(0)  =  0, 

as  seen  by  summing  m(0),  q(0)  and  s(0)  in  Table  I.  This  state  has  the  same  feature 
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of  a  photon-trapping  state  for  the  case  of  two  atoms  in  free  space.  For  the 
initial  state  (e),  each  atom  is  itself  in  an  in-phase  combination  of  excited  and 
ground  atomic  levels.  On  the  average,  each  atom  contributes  a  half  unit  of  excita¬ 
tion;  therefore  there  is  one  total  unit  of  excitation,  m(0)=l,  and  thus  q(0)  =  0. 

The  enhancement  is  not  as  strong  as  in  case  (c),  owing  to  the  fact  that  the  coherent 
combination  is  not  for  the  system  but  for  the  individual  atoms  only,  giving 
s(0)=  1.  By  the  same  token,  the  trapping  effect  of  a  photon  is  small  as  well, 
so  that  s(0)=  -1  for  the  initial  state  (f),  in  which  one  atom  is  in  an  in-phase 


and  the  other  is  in  out-of-phase  combination  of  the  excited  and  ground  atomic 
levels . 


For  the  last  case  (g),  each  atom  is  in  a  random-phase  combination  of  the 

excited  and  the  ground  atomic  levels.  The  excitation  stored  initially  in  the 

atomic  system  is  the  same  as  in  cases  (e)  and  (f),  since  the  radiation  energy 

has  nothing  to  do  with  phase  relation  in  the  atomic  states.  However,  there  is 

no  coherence  resulting  from  the  interaction  with  the  radiation  because  of  a 

null  effect  of  random  phases.  In  summary,  in  comparison  with  the  corresponding 
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cases  of  two  different  atoms  in  free  space,  the  emmission  rates  of  such  an  atomic 

pair  with  dipole  transitions  perpendicular  to  a  metal  surface  are  doubled  in 
magnitude  because  of  the  induced  images  participating  in  phase  with  the  emission 
from  the  source  atoms,  giving  rise  to  a  "superradiati ve"  enhancement  in  the  emis¬ 
sion  rate.  This  effect  is  completely  due  to  the  presence  of  the  metal  surface. 

As  a  final  remark,  we  note  that  the  method  of  solution  presented  in  this 
paper  is  based  on  the  assumption  that  the  interaction  of  electromagnetic  radiation 
and  the  atoms  is  confined  to  a  frequency  uQ  between  and  u^,  with  a  range  roughly 
equal  to  6u).  Therefore,  the  method  is  valid  only  under  the  condition  A  <  5u. 
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To  see  the  time  dependence  of  r-j  ( t)  =  <  R-j  ( t) >  and  ( t )  =  <R2(t)>,one  has  to 
examine  the  equations  of  motion  for  R-j ( t)  and  R^( t) .  A  typical  term  in  R^  and  R^, 
as  given  by  Eqs.(3.39)  and  (3.40),  is  etc. etc.  with  i/ j,  where  i  and  j  represent 

•  *  J  J 

the  double  indices  Ay  or  Iv.  The  equation  of  motion  is 


W  &  <clc1cjcJ>  ■  t4cicjcj-Heff]  - 


-Kg  I  I  [eje,  -  KgMa^cJc^cJcj.Cj] 

•*9“cfo[cjcj-cKr  (A>> 

Making  use  of  the  anti  commutation  and  commutation  relations  (2.1)  and  (2.2),  we  have 


[‘mV'M  *  cfci*ciJ  *  t4'ci]ci  *  -  s 

[ctci ,c|]  =  ct[Ci,ct]  +  [ct,ct]c.  =  ct 


so  that  Eq.(Al)  becomes 

*  3T  (c(cicIcj>  *,|5«a+6(c,c]c>4c1cj) 

Ot  p 

'*9aB(C<CJCJ+Clc1cj)*o6 

4* 

Next,  we  substitute  Eq .(3.16)  and  its  Hermitian  conjugate  for  a^g  and  a^  in 
the  above  equation  and  take  the  expectation  value  to  obtain 

S  -«"*92P(<a0)=c{c1c]cj> 


sWiVjY  =  -Wpf^cfocjc^  -  -ZlKcfocJcj: 
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1.  Two  two-level  atoms  with  different  resonance  frequencies  and 
located  near  a  perfectly  conducting  metal  surface.  Both  electronic 
transition  dipoles  are  assumed  to  be  perpendicular  to  the  surface 
plane,  (a)  direct  photon  transmission;  (b)  photon  transmission  via 
reflection  from  the  surface. 

2.  Normalized  emission  rate  for  the  initial  state  in  Eq .(4.1). 

3.  Normalized  emission  rate  for  the  initial  state  in  Eq .(4.2) 

4.  Normalized  emission  rate  for  the  initial  state  in  Eq.(4.3). 

5.  Normalized  emission  rate  for  the  initial  state  in  Eq.(4.4). 

6.  Normalized  emission  rate  for  the  initial  state  in  Eq . (4.5) . 

7.  Normalized  emission  rate  for  the  initial  state  in  Eq .(4.6). 

8.  Normalized  emission  rate  for  the  initial  state  in  Eq .(4.7). 
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